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5 : Abstract 

^' The fermion mass spectrum is studied in the quenched approximation in the strong 

^ . couphng vortex phase (VXS) of a globally U (1) l<Si^ (l) r symmetric scalar-fermion model 

in two dimensions. In this phase fermion doublers can be completely removed from the 

'^>( I physical spectrum by means of a strong Wilson- Yukawa coupling. The lowest lying fermion 

Vh I spectrum in this phase consists most probably only of a massive Dirac fermion which has 

charge zero with respect to the U{1)l group. We give evidence that the fermion which is 

charged with respect to that subgroup is absent in the VXS phase. When the U{1)l gauge 

fields are turned on, the neutral fermion may couple chirally to the massive vector boson 

state in the confinement phase. The outcome is very similar to our findings in the strong 

coupling symmetric phase (PMS) of fermion-Higgs models with Wilson- Yukawa coupling in 

four dimensions, with the exception that in four dimensions the neutral fermion does most 

probably decouple from the bosonic bound states. 
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1 Introduction 

A non-perturbative formulation of a chiral gauge theory on the lattice has proved to be a 
difficult issue. In a chiral gauge theory naively transcribed to the lattice, each fermion is 
accompanied by fifteen 'doubler fermions' where eight of these couple as mirror fermions and 
spoil the chiral couplings. One way to deal with this problem is to decouple the unwanted 
extra species by rendering them very heavy. 

The standard Wilson mass term, which is known to remove the doublers in the case 
of vector-like theories on the lattice, obviously breaks gauge invariance of the chiral gauge 
theory. A proposal to overcome this difficulty is the so-called Wilson- Yukawa approach [0, ||] 
which has recently received a lot of attention. Instead of a standard Wilson mass term one 
uses a so-called Wilson- Yukawa term which contains the Higgs fields in a way that it is 
manifestly invariant under the chiral gauge transformation. The goal of prime importance is 
now to try to decouple the unwanted species doublers by means of a strong Wilson- Yukawa 
coupling and to give them a mass of the order of the cutoff. It is well known that such a 
decoupling is a non-trivial and non-perturbative issue. For recent overview articles on the 
Wilson- Yukawa approach see refs. [3-9]. 

Recently extensive investigations of globally G^^C^Gr symmetric (with Gl^h =SU(2), 
U(l)) fermion- Higgs models with Wilson- Yukawa coupling in four space-time dimensions 
have shown that it is rather unlikely that this method leads to the desired non-perturbative 
formulation of the standard model on the lattice. The reasons for this may be summarized 
as follows: The phase diagram contains, apart from weak coupling symmetric (PMW) and 
broken (FM(W)) phases where fermion masses behave according to perturbation theory, also 
strong coupling symmetric (PMS) and broken (FM(S)) phases where these masses exhibit a 
non-perturbative behavior. The Wilson- Yukawa approach fails in the weak coupling phases 
because there the masses of the doubler fermions are restricted by the triviality of Yukawa 
couplings and cannot be made sufficiently heavy ||10|. They remain as additional particles 
in the spectrum. 

On the other hand, in the strong coupling phases PMS and FM(S) the doubler fermion 
states can be removed completely from the particle spectrum by making them as heavy as 
the cut-off. In contrast to the weak coupling region, fermions become massive also in the 
PMS phase, a situation which differs already from that in the fermion-Higgs sector of the 
perturbative standard model where the fermion mass vanishes in the symmetric phase. Sym- 
metry considerations show that the particle spectrum in the PMS phase can a priori contain 
a fermion which is neutral with respect to the G^ group and a fermion which is charged with 
respect to that group [11-19]. The existence of the neutral fermion was confirmed by the 



good agreement of analytic and numerical calculations |]T3, O, [T^. On the other hand a 1/w 



expansion {w is the Wilson- Yukawa coupling) of the charged fermion propagator [2C] and a 



numerical investigation of the fermion spectrum [Q gave evidence that the charged fermion 



does not exist as a particle in the spectrum, though this is not generally accepted ||2T| . Under 
the assumption that the charged fermion is absent one can show that the coupling of the 
neutral fermion to Higgs and gauge fields vanishes as a power of the lattice spacing a and 



the neutral fermion becomes non-interacting in the scaling region |13|. It was also argued 



that the renormalized Yukawa coupling in the FM(S) phase vanishes as a power of a rather 



than logarithm of a [0, [14 



In this paper we extend the investigations to a U(1)l(S'U(1)r symmetric fermion-Higgs 



model with Wilson- Yukawa coupling in two dimensions. The important advantage of an 
investigation in two dimensions is that the simulations can be carried out on lattices of large 
linear dimension, enabling one to achieve large correlation lengths for the scalar fields with 
better control over finite size effects. The numerical data are of much superior quality than 
in the four dimensional case. 

In two dimensions spontaneous breakdown of a continuous symmetry cannot occur be- 
cause of the Mermin- Wagner- Coleman theorem P^. In spite of that, fermions are observed 



to acquire a mass also in two dimensional fermion-scalar models. For example a 1/A^ ex- 
pansion in the Gross Neveu models with a continuous chiral symmetry (which can be also 
viewed as fermion-scalar models after the introduction of an auxiliary scalar field) shows that 
the fermion mass does not vanish |^|. At a first glance this appears to be contradictory. 



The contradiction could be resolved by expressing the action in terms of new fermionic fields 



which can have a mass term without destroying the original chiral symmetry pj]. Inter- 
estingly, the existence of massive fermions in the strong coupling symmetric phase PMS in 
the four dimensional models with Wilson- Yukawa coupling may be viewed in a similar way 



n, O. The new fermionic variables are here the above mentioned neutral and charged 



fermion fields which transform vectorially under the original chiral symmetry group and al- 
low therefore for the construction of invariant mass terms. 

The phase structure of the two dimensional model in the quenched approximation is 
similar to the one found before in the four dimensional models. The analogues of the weak 
coupling phases FM(W) and PMW are respectively a weak coupling spin-wave SW(W) and 
vortex VXW phase. The strong coupling phases FM(S) and PMS get replaced by strong 
coupling spin-wave SW(S) and vortex VXS phases. We find that fermions are massive in 
the strong coupling phases VXS and SW(S). The existence of the strong coupling phases is 
recently confirmed also by an investigation of a two dimensional U(l) fermion-Higgs model 
with dynamical naive fermions. There are, however, some indications favoring the absence 
of the VXW phase ^5\. The weak coupling spin-wave SW(W) would in this case extend 



down to zero Yukawa coupling. If this is correct, the VXW phase has to be regarded as an 
artefact of the quenched approximation. 

In this paper we study the fermion spectrum in the quenched approximation in the strong 
coupling vortex phase VXS whose existence is guaranteed in the full model with dynamical 
fermions. In this phase the species doublers can be completely removed from the spectrum 
and the generation of the fermion masses is based on the same mechanism as in the PMS 
phase of the four dimensional models. Also here the fermion spectrum can a priori consist 
of a neutral and a charged fermion. Similar to the four dimensional models we find strong 
indications for the absence of the charged fermion. The spectrum consists then only of the 
neutral fermion and the scalar particles. We shall show in sect. 7 that the neutral fermion 
may exhibit a chiral coupling to the vector boson state in the confinement phase after the 
gauge interactions are turned on. 

The outline of the paper is as follows: In sect. 2 we introduce the model and describe its 
phase structure in the quenched approximation. Sect. 3 deals with the fermion spectrum in 
the strong coupling phases SW(S) and VXS. In sect. 4 we present the results of a leading 
order hopping expansion for the neutral and charged fermion propagators. After giving a 
brief report on the technical details of the numerical simulations in sect. 5, we compare in 
sect. 6 the numerical results for the rest energies of the charged and neutral fermion with 



those obtained from the hopping expansion. Based on the results of sect. 6 we discuss in 
sect. 7 different outcomes for the physics in the VXS phase. A brief conclusion is given in 

sect. 8. 



2 The model and its phase diagram 

The model of interest is given by the following gauge invariant euclidean lattice action in 
two dimensions: 

S = ^^ , ^ = >Cgauge + >Cgcalar + ^fermion (2- 1) 



with 
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+ $:*x + A($:$x - 1)' , (2.3) 
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where in the second line of eq. ( p.4|) we have included the Wilson- Yukawa coupling with 

strength w. Besides the Wilson- Yukawa term we also included a usual Yukawa term. The 

symbols D 

defined by the relations D'^'^ 

and d~ 



d^ and d denote the covariant and normal lattice derivatives which are 

U* -vTf - B+ = n^ 



Ufj,X^X+jl 



^x, D-^, 



^x 



\U=l 



D-\u- 



=1. The symbols 7^, /x = 1,2 denote the two dimensional 7 matrices and 
^R,L = \{^ =1=75) with 75 = — i7i72 are the right and left-handed chiral projectors. The field 
Ui2x = UixU2x+iU*^,^U2x is the usual plaquette variable in the Wilson action, g is the gauge 
coupling, K is the hopping parameter for the scalar field and A the quartic coupling. If we 
succeed in removing the species doublers from the particle spectrum by means of a strong 
Wilson- Yukawa coupling w, one would expect that the lattice lagrangian defined by eq. ( p]) 
reproduces in the continuum limit the target model given by the lagrangian 



-^0 — -7Ffj_^Ff^,y + Dfj_(f 



'D''(j) + ml(t)*(t) + Ao 



+ij{It)PL + ^Pr)^^ + Z/oV^(0Pi? + 0*Pl)V' 



(2.5) 



where F^y{x) 



d^Aj 



X] 



dyA^{x) and D^ = d^ — igoA^{x). We will demonstrate in this 
paper that this naive expectation is not correct and that the effective lagrangian which 
describes the physics in the VXS phase differs substantially from eq. ( p.5| ). The continuum 
fields ip{x), (j){x), Af^{x) in eq. (^75|) are related to the corresponding lattice fields in ( |2.1|) by 
the transformations 
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-^(j){x) , U/^x = exp{-iagoA 
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The coupling parameters in eq. 
the relations 

o l-2A-4/t 



can be expressed in terms of the lattice couplings by 

A 
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We note that in the continuum formulation the Yukawa couphng, the gauge couphng and the 
quartic couphng carry a dimension, whereas the gauge and scalar fields are dimensionless. 
Throughout this paper we will study the model in the limit A ^ oo which implies that the 
scalar fields become radially frozen, i. e. $*$x = 1- 



The lattice lagrangian in eq. ( |2.1| ) is invariant under the local gauge transformations of 
the form '^l,x ^ ^l,x'^l,x, "^l.x ^ "^l.^^I,^, ^x ^ ^l,x^x, U^^ -^ VLl,xU^:,VlI^^^^, with 
^L,x £U(l)j;,. The model is furthermore also invariant under the global gauge transformations 
^R,x -^ ^R^R,x, '^R,x -^ '^R,x^*R and $^ ^ <l>^fi|j with Qr gU(1)r, At y = the action 
( p.l|) possesses a shift symmetry for the right-handed fermion fields, 

'^R,x -^ ^R,x + eR , '^R,x -^ '^R,x + Ir , (2.8) 

where eR and Ir are the constant shifts in the right-handed Weyl spinors. This symmetry 
guarantees that the mass of the fermion with the quantum numbers of the "^r fermion field 
vanishes at y = for all values of n and g 



Although we have included for later convenience in ( [^.ID also the interactions to the U(1)l 
gauge fields, in our numerical work, however, we will restrict ourselves to the global limit 
g = where Uf^x = 1 and the local U{1)l gauge symmetry turns into a global one. In this 
case eq. ( p.3| ) (with A = oo) reduces to the lagrangian of the XY model in two dimensions. 



We shall furthermore study the model (|2.1| ) in the quenched approximation where the effects 
of the fermion determinant are neglected. The use of the quenched approximation for the 
investigation of the fermion spectrum in strong coupling phase VXS is justified since this 
phase was established also in the full model with dynamical fermions 



Next we turn to the phase structure of the model in the quenched approximation. The 
XY model is known to have a phase transition at k = Kc ~ 0.56 which separates a vortex 
(VX) phase {k < Hc) with finite scalar correlation length from a spin-wave (SW) phase where 
the scalar correlation length is infinite p^. We note that the spectrum in the VX phase 



consists of two scalar particles which have the same mass. In the quenched approximation 
fermions do not have a feedback on the scalar sector and Hc is independent of y and w. 
A phase transition may, however, occur in the fermionic sector. Such a phase transition 
or 'crossover' was discovered before in four dimensional models with Wilson- Yukawa cou- 
pling at y + Aw ^ \f2 |2^, ^. It separates the weak coupling phases PMW and FM(W) 



(y + 4w ~ \pi) from the strong coupling phases PMS and FM(S) (y + 4w ~ \pi)- The 
existence of this 'crossover' has manifested itself in a different behavior of the fermion mass 
as a function of n in the weak and strong coupling regimes. 

In order to monitor the fermionic phase structure of the quenched model in two dimen- 
sions we have computed the mass mp of the ^r fermion by fitting the (^ r^ rj propagator 
to a free Wilson fermion propagator ansatz (for more information about the technical details 
see sect. 5). In fig. 1 we have displayed this mass as a function of k for several values of the 
Yukawa coupling y and for the special case of naive fermions (w = 0). As in four dimensions 
the fermion mass shows a qualitatively different behavior at small and large values of y. It 
decreases when approaching the SW-VX phase transition in the SW phase for y < y* ^ 1. 
On the contrary the mass is observed to increase when the value of k, is lowered for y > y*. 
This increase is seen to continue into the VX phase and fermions become massive in that 
phase. A similar behavior of the physical and doubler fermion mass as a function of k is 
observed also for w > 0. From the different behavior of the fermion mass we can localize 



Figure 1: The mass of the "^r fermion as a function of y for several values of k at w = 0. 
The position of the phase transition in the thermodynamic limit is given by Hc ~ 0.56. The 
vertical arrow below the abscissa indicates the position of the peak in the susceptibility on the 
48 X 48 lattice which lies clearly below Kc- 

the position of a 'crossover'. This position is within the precision of our resolution indepen- 
dent of K and is approximately given by the relation y + 2w = 1. As in four dimensions 
the 'crossover' sheet splits the VX and SW phases into strong (S) and weak (W) coupling 
regions, which we will denote in the following by VXS, SW(S) {y + 2w ~ 1) and VXW, 
SW(W) (?/ + 2w ~ 1). As we mentioned already in the introduction, the VXW phase does 
not seem to be present in the full model with dynamical fermions and could be an artefact 



of the quenched approximation [^. However, this does not affect our investigations in the 
VXS phase. 

Our numerical results on the 48^ lattice show that the fermion mass stays non-zero ev- 
erywhere in the SW phase and also in the VXS phase though the chiral symmetry cannot 
be broken according to the Mermin- Wagner- Coleman theorem. Studies on different lattices 
show that the finite size dependence of the fermion mass is extremely small and it appears 
very unlikely that the fermion mass could vanish in the infinite volume limit. The existence 
of a massive neutral fermion in the strong coupling phases SW(S) and VXS is indeed strongly 
supported by the good agreement between analytic calculations which are based on strong 
coupling expansions and the results of the numerical simulation. We will report more on 
these results in the following sections of this paper. 



The fact that fermions may be massive within the SW phase although the chiral symmetry 
is unbroken in that phase was to our knowledge first explained in ref. [^]. The basic idea is 



as follows: The original action in terms of the \E' fields does not provide a correct description 
of the physics in the SW phase. It can, however, be rewritten in terms of new fermionic field 
variables, which transform vectorially under the original symmetry transformations thereby 
allowing for the construction of a chirally invariant mass term. This new form of the action 
may give a more appropriate description of the physics in this phase (in the sense of a weak 
coupling expansion) if the fermions are indeed observed to be massive. 

3 The fermion spectrum in the strong couphng phases 

Since the \J{1)l®^{1)r symmetry is unbroken everywhere, the states excited by the fields 
"^L and \E'ij need not be the same since the corresponding fields carry different quantum 
numbers under the unbroken symmetry group. We will refer to the \I/r = \E'j^ field as the 
neutral (n) fermion field since it has charge zero under the U(1)l group {qi = 0). The 
"^L = ^L field will be called charged (c) fermion field since it has charge one under the local 
U(1)l gauge group (g^ = 1). When discussing the phase diagram in the previous section we 
have already mentioned that the numerical results give strong evidence that the \l/ij fermion 
becomes massive both in the SW(S) and VXS phases. This will be substantiated later by the 
good agreement of the numerical results with a hopping expansion. In order to describe this 
situation the "^^ field may be regarded as the right-handed component of a massive Dirac 
fermion field \Ef("). A possible choice for "^^^ is the composite field ^*"^l which transforms 
in the same manner as "^ji- The neutral Dirac fermion field may then be written as 

^(n) = ($*p^ + p^)^ ^ ^(") = "^($p^ + p^) . (3.1) 

Along the same lines we may also introduce a charged Dirac fermion field 

^W = (p^ + $P^)^ , "^^'^ = "^(Pr + $*Pi) . (3.2) 

The fields '^^'^^ and \E'*^"^ transform vectorially under U(l)i and U(l)ij respectively. 

On a finite lattice the long range fluctuations in the SW phase cause a non-vanishing 
value of the magnetization M which may be defined by the relation 

^=(^E'^=^)lrot, (3-3) 

where V is the lattice volume. The index "rot" means that each configuration is rotated 
such that ^ J2x ^x points into a given direction. This rotation is necessary since on a finite 
lattice the magnetization M is drifting through the group space and when averaging over 
many configurations one would get zero. A measure for the magnetization in the SW phase 
may then be defined by f = |M|. Since there is no spontaneous symmetry breaking in two 
dimensional systems, this quantity has to vanish in the infinite volume limit V ^ oo. Some 
typical values for w on a 48 x 48 lattice are given by v{k) = 0.3545(44), 0.5565(22), 0.6734(12) 
for K = 0.48, 0.52, 0.60. Even on large lattices (e. g. 400^) v is clearly non-zero in the SW 
phase and increases when raising the values of k. Therefore, on a finite lattice the situation 
in the SW phase is very similar to the broken phase in the four dimensional model where the 
U(l)L®U(l)ij symmetry is broken to the diagonal subgroup \J{1)l=r- As a consequence the 



fields \1/''"'^ and '^^^^ appear to behave almost as equivalent interpolating fields. Indeed, the 
numerically found rest energies obtained from the neutral and charged fermion propagators 
coincide in that phase within the statistical errors. In the infinite volume limit, however, the 
two rest energies are expected in general to be different for w > 0. 

4 Hopping expansion for the neutral and charged 
fermion propagators 

An appropriate method to try to find analytic expressions for the neutral and charged fermion 
propagators in the strong coupling phases is the hopping expansion. The hopping expansion 
deals only with the fermionic integration in the path integral, the bosonic integration has 



to be performed e. g. by a 1/d expansion |jT5| or by numerical simulations. When starting 



from the lagrangian ( p^.4| ) with the single-site Yukawa-coupling and expanding the Boltzmann 
factor in the path integral in powers of the hopping parameter a = l/{y + 2w) one comes 
upon cancellations of the type $*$ — > 1, emerging from the single-site terms and the one-link 
terms. The hopping expansion becomes more transparent after removing the $ fields from 
the single-site Yukawa term by means of a unitary transformation to the one link terms. Two 



transformations of this type are given by the inverses of eqs. (|3.1|) and (|3.2| ) which express the 
neutral and charged fermion fields in terms of the original \E' fields. The associated jacobian 
for these transformations is in both cases equal to one since the scalar field is radially frozen. 
Replacing \E' and "* in (|]4D by \E'(") and "^^""^ gives 



^F = W 



2.^1 



wSr)$*)7,p: + D-)i<^¥^^) +¥^\{d:^ + d-)¥-^ 



+^^(")vi/W _ !f!^(") J2 ^,t^;^^"^ • (4-1) 

This substitution transforms the Yukawa term into a bare mass and the Wilson- Yukawa 
term into a free Wilson term. Expressing the ^ and ^ fields in terms of the "^^^^ and ^ 
fields leads to 



1 



2 



7^., rr>+ I n-^1Tf(c) , ^Ti7W^^„, /a+ , ^-\r^*,TM)^ 



2 

+2/^(^)^W _ -(^^"^c^) ^ d+d-{<l>*¥'^) . (4.2) 

The lagrangian ( [4.1| ) has a shift symmetry ( ^.81 ) in terms of the neutral field because 
\E'^ = \E'/j. Such a symmetry, however, is absent for the action ( [4.2| ) in terms of the "^^^^ 
fields. This different behavior under the shift symmetry holds also if the local 11(1)^ gauge 
interactions are turned off. It makes therefore sense to distinguish between the charged and 
the neutral fermion also in this globally symmetric U(1)l(8>U(1)/j theory. 



Using the lagrangians ( [4.1| ) and ( [4.2| ) one finds to lowest order in a the following expres- 
sions for the charged and neutral fermion propagators in momentum space. 



2 2 

[z-^Pl + PbMv + 2w-wY^ cos k^)z-^ + ^ E 7/. sin k^]-\z-^PR + Pl) 



(4.3) 



and 



S^''\k) = /l^^W^^'^^e'^^""'^) 



2^,y 
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(Pi + z-^Pi?)[(z/ + 2w)z-^ - wzJ2 cos k^ + I J2 l^.smk^]~\PR + z'^Pl) 



(4.4) 



Here 



z^ = (Re($:f/^.<l>,.+^)) (4.5) 

is the scalar field link expectation value. This quantity has a non-vanishing value in the VX 
as well as SW phase. 

From the expressions (|4.3|) and (|4.4|) one can read off expressions for the fermion masses. 
For the masses of the neutral fermion and its species doublers we obtain 

mP ^ yz-^ , mS^^ ^ mj^^ + 2wlz~^ , / = 1, 2 , (4.6) 

where / is the number of momentum components equal to vr in the two dimensional Brillouin 
zone. From eq. ( |4.(j| ) we can read off an effective Wilson r-parameter: 

r(") ^ wz-^ . (4.7) 

As we shall see later, these expressions are in good agreement with the numerical simula- 
tions in the VXS phase. This was also found to be the case in the four dimensional models 
|ll,[T7|, HI. 



In agreement with the shift symmetry [^ the mass nip of the physical fermion is seen 
to vanish in the limit y ^ 0. The doubler fermions, however, are non-zero within the strong 
coupling phases, since 2wlz^^ ^ everywhere in this region. This means that in the contin- 
uum limit the species doublers for the neutral fermion decouple from the particle spectrum 
within this phase. 



Similar formulas can be obtained from eq. (|4.4| ) for the masses of the charged fermion 
(assuming it exists for the moment) and its species doublers 

rrip ~ (y + 4:w)z~'^ — Awz , m^^ ^ nip + 2wzl , / = 1, 2 . (4.8) 

The effective Wilson parameter is now given by 

r^") ^ wz . (4.9) 

The discussion in the following section will show that the formulas ( ^4.8| ) and ( |4.9| ) for the 
charged fermion are in disagreement with the numerical results. A calculation of higher 
order terms in four dimensional models showed that they appear to be small for the neutral 
propagator S^"'\ but not for the charged propagator S**^^^ ||18|, O. On the basis of our 



numerical results we will give in sect. 7 an argument why the hopping expansion to lowest 
orders leads to a wrong result for charged fermion propagator. In sect. 6 we will compare 
the analytic formulas that have been derived in this section with the results of the numerical 
simulation. 



5 Details of the numerical simulation 

The neutral and charged fermion propagators were determined by inverting the correspond- 
ing fermion matrices on a set of uncorrelated scalar field configurations which were generated 
by means of the refiection cluster algorithm [^ for the XY model. We have computed the 
fermion propagators in coordinate space. As an example we give the expression for the RR 
component of neutral fermion propagator 

41(i) = (^ E '^^Ln^tL/''^''''''^) . i=\^2- y2\ , (5.1) 

where t = 1, . . . , T. The symbols L and T denote here and in the following the spatial and 
time extent of a rectangular lattice. The physical fermion propagator is obtained for pi = 
and the propagator of the lowest lying doubler fermion for pi = ir. The fermion and the 
doubler fermion propagators were computed for all four L-R combinations. We have used 
for the fermion fields periodic boundary conditions in the spatial direction and anti-periodic 
boundary conditions in the time direction. The scalar fields had periodic boundary condi- 
tions in all directions. 

For the neutral fermion propagator we have inverted the fermion matrix on typically 
1000-2000 scalar field configurations. A problem which we were confronted with in the four 
dimensional models was the large number of matrix inversions which was required to obtain 
a stable signal for the charged fermion propagator in the PMS phase. In the two dimensional 
models it is possible to enlarge the statistics in the VXS phase on relatively large lattices 
(e. g. 32^) by an order of magnitude. For the determination of the charged fermion prop- 
agator we have inverted the fermion matrix on typically 1000-5000 and deeper in the VXS 
phase on 20000 scalar field configurations. 

Most of the results were obtained on a 32 x 32 lattice. In order to estimate the finite size 
effects we varied L at a particular point in the VXS phase from 16 to 64 while keeping T 
fixed at 64. 

We have fitted the neutral and charged propagators at zero spatial momentum to the 
free Wilson fermion ansatz 



2\/l + 2rimi + mf 



exp{-Eit) + Cexp{-Ei{T-t)) _ , expj-E'it) + Cexp{-El{T - t)) 

l + exp{-EiT) ^^ ' l + exY>{-E[T) 



(5.2) 



where C, = 1(— 1) for the RR and LL (RL and LR) components. This ansatz holds only for 
ri < 1, for r/ > 1 the oscillating factor (—1)* has to be omitted. 

From this fit we obtain the numerical values for the rest energies Ei and E'l of the fermion 
and its 'time doubler' and for the wave-function renormalization constant Z. The subscript 
/ is for the physical fermion propagator (pi = in eq. ( |5.1| )) and 1 for the doubler fermion 



propagator {pi = vr). The masses mi and the Wilson parameters ri for which we obtained 
expressions in the previous section are related to the rest energies Ei and E'^ by the lattice 



Figure 2: The physical charged fermion propagator Sfil{t) (pi = 0) as a function of the 
time coordinate t for (k, y, w) = (0.45, 0.3, 2.0). The lattice is size is 32 x 64. The solid line 
was obtained by fitting Sj!il{t) to the free Wilson fermion ansatz in eq. (S.i). 



dispersion relations 



yl + 2rimi + mf + ri + mi Jl + 2rimi + mf + ri + mi 

e ' = ^ -— , e ' = ^ . (5.3) 

1 + n l-ri 

The effective Wilson parameter r/ could be in principle a function of / = 1,2. The numerical 
results, however, show that r/ is independent of / (therefore we will subsequently use the 
notation r = vi), which supports the interpretation of the numerical results in terms of the 
free fermion formula. 

We find that the rest energies obtained from the four chiral components Sj^l{t), 5')^(t), 
S^l{t) and Si^^{t) for pi = 0, vr agreed always within the statistical errors. The same holds 
also for the rest energies determined from the four chiral components of the charged fermion 
propagator. In the following we will use the notation Ep and Ep for the rest energies of 
the neutral and charged physical fermion and similarly Ej^ and E^ for the rest energies of 
the lowest lying doubler fermions. 

In fig. 2 we have displayed as an example the charged fermion propagator S^l{t) for 
Pi = as a function of t. In this example we have chosen L = 32 and T = 64. The solid 
curve was obtained by fitting the numerical data to the free Wilson fermion ansatz ( |5.2| ). 
The high quality of the numerical results for the propagators allowed for an accurate deter- 
mination of the rest energies E^ , E^ , E^ and E^^\ 



For the considerations in sect. 7 we have to know also the numerical values for the 
rest energy £"$ of the scalar particles in the VX phase. This rest energy was determined 
numerically by fitting the scalar propagator in momentum space defined by 

GM = l^^Y.^:%e^'^^-y^^ (5.4) 

to a free scalar propagator ansatz 

where p^ = 2X]^=i(l — cosp^) is a lattice equivalent of the momentum squared in the 
continuum. The rest energy ii^$ of the scalar particles in the VX phase is obtained from the 
lattice dispersion relation at zero momentum, coshE'^ = 1 + m|/2. 

6 Comparison of the numerical results with the hop- 
ping expansion 

In this section we are going to compare the numerically found values for the rest energies 
Ep and Ep of the neutral and the charged fermion and the corresponding rest energies 
E^^ and E^/ for the lowest lying doubler fermions with the analytic expressions which result 
from the formula ( p. 3D after inserting eqs. ( |4.6|) , ( |4.7| ), ( [4.8|) and ( |4.9| ) from the lowest order 
hopping expansion. For z^ we use the numerical value measured on the same lattice which 
we are using for the determination of the fermion propagators. In fig. 3 we have displayed the 
numerical values for Ep\ Ej^ , Ep and Ep^ as a function of y for k = 0.4 and w = 2.0. The 
coupling parameter values we have chosen here lie well inside the VXS phase. The results 
from the hopping expansion are represented by the curves. The dashed, solid, dash-dotted 
and dotted lines correspond respectively to the rest energies Ep , E}^, Ep and E^) . The 
figure shows that the agreement between the numerical result and the analytic prediction is 
quite impressive for the rest energies Ep and E^ while the curves for Ep and Ep) exhibit 

(c) 

a strong deviation from the numerical results. In the case of Ep the deviation is larger than 

(n) 

a factor two. The figure shows furthermore that Ep appears to vanish in the limit y ^ 0, in 
agreement with the shift symmetry mentioned before, whereas Ej^ stays non-zero for all dif- 
ferent values of y which implies the decoupling of the species doublers of the neutral fermion 
in the continuum limit. Also the numerical values for Ei^ are larger than one for all values 
of y with no indication of dropping in the limit y ^ 0. Thereby also the species doublers 
of the charged fermion can be completely removed from the physical spectrum. Provided 
it exists at all as a particle in the spectrum, the charged fermion is massive in the VXS phase. 

(n) (c) 

In fig. 4 the rest energies Ep and Ep are plotted as a function of k where the coupling 
parameters y and w were fixed to 0.4 and 2.0. The analytic results for Ep and Ep are 
represented also in this figure by the dashed and dash-dotted lines. We find again that 
analytic expressions from the hopping expansion provide a good description of the rest 
energy of the neutral, but not of the charged fermion. The other details in this figure will 
be explained in the next section where we will develop, on the basis of the results of this 
section, two different scenarios for the physics in the strong coupling region. 



Figure 3: The rest energies Ep , Eq , Ep and E^ as a function of y for k = 0.4 and 
w = 2.0. The dashed, solid, dash-dotted and dotted curves represent respectively the analytic 
results for Ep , E}^ , Ep and Ej^ obtained from the hopping expansion. 

7 Scenarios for an effective field theory in the strong 
couphng regime 

The results of the previous section showed that the analytic results deduced from the la- 
grangian ( [4.1|) are in good agreement with the numerical results for the rest energies Ep 
and Ej^ , whereas the lagrangian ([4.2|) led to expressions which were in a strong disagree- 



ment with the numerical data. This suggests that the physics in the strong coupling region 
is well described by the lagrangian ([4.1| ) in terms of the neutral fermion fields. The charged 
fermion fields '^^'^' = ^'^^'^' and \l/ = v]/ $* can then be regarded as composite fields and 
the charged fermion, provided it exists at all in the particle spectrum, can be considered 
as a bound state of the scalar particle and the neutral fermion. The question arises now 
whether the interactions in eq. i\i.l\) can produce such a ^-vp^") bound state. In the four 
dimensional model the scalar fields $ carry a dimension of a mass and as a consequence of 
this the four-point coupling 

lj:i^P^*h,{D; + D;){m'£^) (7.1) 

with d = A vanishes in the classical continuum limit like a^ which makes the formation of a 
$_^(") bound state already very unlikely. Indeed a 1/w expansion [^ and the numerical 



Figure 4: The rest energies Ep and Ep plotted against k for y = 0.4 and w = 2.0. The 
dashed and dash-dotted curves represent the analytic results from the hopping expansion for 
Ep and Ep . The solid line gives the result for the sum Ep + £"$. The error bars for 
Ep + Eq, are much smaller than the symbol sizes for Ep . 



simulations |T^ gave strong evidence for the absence of the charged ferniion in the particle 
spectrum of the PMS phase. In the case of the two dimensional model the scalar fields are 
dimensionless and for d = 2 the coupling ( [7. ID does not vanish as a power of the lattice 
spacing a. Therefore the formation of a $-\['(") bound state appears at the first glance to be 



more favored than in the four dimensional model. If the four-point interaction (|4.1| ) is strong 
enough to produce a $-\[f(") bound state we expect to find the following relation among the 
rest energies of the neutral and the charged fermion and of the scalar particles in the VXS 
phase 

EP = EP + E^ + eB, eB<0, (7.2) 

where the quantity e^ denotes the binding energy of the <|)-\[f("') bound state. This relation 
implies that the charged fermion could only exist as a particle at a point in the coupling 
parameter space where the rest energies Ep , Ep and es scale simultaneously to zero. This 
can happen only at the point k = Kc, y = 0, since only there Ep and Eq, can vanish simul- 
taneously. 

Scenario A: Let us assume now for the moment that e^ scales to zero like Ep and Eq, in 
the limit k y Kc, y ^ and the charged fermion exists together with the neutral fermion as 
a Dirac fermion in the particle spectrum. The coupling of the '^^ field to the gauge fields is 
necessarily vectorial because the charged fermion is massive (m^'^^ > 0). The model we end 



L 


hi^ 


eP 


E^^ 


es 


16 


0.126(3) 


0.275 


0.368(13) 


-0.033(16) 


32 


0.115(3) 


0.277 


0.374(9) 


-0.018(12) 


48 


0.120(5) 


0.276 


0.379(12) 


-0.017(17) 


64 


0.119(5) 


0.277 


0.389(11) 


-0.007(16) 



Table 1: The rest energies E^, Ep , Ep and es as a function of the spatial extent L of a 
lattice with volume L x 64 obtained at the point {K,y,w) = (0.45,0.3,2.0). The error bars 
for Ep were omitted since they are smaller than 0.001. 



up with in the strong couphng VX phase is clearly different from the original target model 
given in eq. (|2.5|), although we succeeded in removing the species doublers from the spectrum. 
One possible form of an effective action in the VXS phase is given by the expression 



CJ' 



v^(") 






'^1 



(7.3) 



where we used the concise continuum notation of eq. ( P~5D and all fields and coupling pa- 
rameters are considered to be effective. The symbol y^ denotes the renormalized Yukawa 
coupling. The failure of the Wilson- Yukawa approach in giving a chiral gauge theory is re- 
lated to the fact that the original fermion fields \l/ij and \l/i combine with composite 'mirror' 
fermion fields xr = 'Pi^R and xl = 0Vl and form the two massive Dirac fields \1/*^"'' and '^^'^^ 
which would couple vectorially to "external" gauge fields. The model (ITTsI) is a special case 
of the mirror fermion model |^T|, when transcribed to the case of two dimensions. It has, 
however, less flexibility in tuning coupling parameters. 

Eq. ( |7.3|) is, however, not the only possible form of an effective action in the strong coupling 
phase. For example it allowed by the symmetries to add a term, which couples the neutral 
fermion chirally to the massive vector bosons in the confinement phase. We will come to 
this in the last part of this section. 



Scenario B: In order to figure out whether the above requirements for the binding energy 
are fulfilled we have determined e^ numerically in a wide range of the bare parameters in 
the VXS phase. The details about the numerical determination of the rest energy E^p were 
given in sect. 5. In the figs. 4 and 5 we compare the rest energy Ep with the sum Ep + Eq, 
which in these graphs is represented by the solid lines. For Ep we have used the results 
from the hopping expansion which are in perfect agreement with the actual data, as we 
reported in sect. 6. The error bars for the sum Ep + E^ are always much smaller than the 
symbol sizes for Ep . Both plots indicate that e^ = for all values of k and y in the VXS 
phase. This suggests that the formation of a bound state is very unlikely. Fig. 5 shows that 
also for the lowest lying doubler fermion the numerical results for the rest energy E]^ are 
nicely represented by the sum Ef^' + Eq, (upper solid line). In order to get an impression 
about the finite size dependence of the energies Ep\ Ep , Eq, and e^ we have computed 
the S^'^\ S^^^ and Gq, propagators at a fixed point in VXS phase on a sequence of lattices 
with size L x 64, where spatial extent L was varied in a range from 16 to 64. The results 
for Ep , Ep , Eq and e^ are summarized in table 1. It can be seen also here that the 
binding energies are very small and even on the smallest lattices almost compatible with 



(c) (c) 

Figure 5: The rest energies Ep and E)-, as a function of y for k = 0.4 and w = 2.0. The 
rest energies are compared with the sums Ep + Eq, and E]^ + Eq, which are represented by 
the lines. The error bars of Ep + Eij, and E}^ + E,^ are much smaller than the symbol sizes 
forE^. 

zero within the quoted error bars. Furthermore one recognizes a systematic trend of le^l to 
decrease when the spatial extent of the lattice is enlarged. These results strongly suggest 
that the signal which we detected in the charged fermion propagators is simply caused by 
a two particle state of the neutral fermion and the scalar particle. We now can understand 
also the fact why the hopping expansion to lowest order for the charged fermion propagator 
is very misleading. To reproduce the inverse propagator of a two particle state an infinite 
number of hopping terms would be needed. Only on the basis of our numerical results we 
can, of course, not completely exclude the possibility that in principle a very small and 
no n- vanishing binding energy might be left over in the continuum limit. We could also not 
find a good field theoretical argument which would rule out the existence of a bound state 
with zero binding energy. Both cases would bring us back to scenario A which we drew up 
first in this section. 

Let us now proceed under the assumption that the charged fermion does not exist as a 
particle in the spectrum. In order to find an effective lagrangian for the model in the VXS 
phase we first rewrite the lattice lagrangian ( [4.1|) in the following way 

1 2 

^F = IE 



V^ L,xl^i^ ^lx^ L,x+fi ^ L,x+il'yiJ.'-^fix ^ L,x) ' l^_R,a:7/i^_R,x+/i ^_R,x+/i7Ai^i?,a;J 






(7.4) 



where we introduced the effective gauge field 



^;, ^ Ku,.^cc+f. . 



(7.5) 



Although the interaction in the first term of eq. ( |7.4| ) appears to be too weak for a formation 
of a <|)-\[f(") bound state, the form ( |7.4|) leaves still the possibility of a chiral coupling be- 
tween the neutral fermion and the effective gauge field U' since this field has dimension one 
according to a naive power counting analysis. This outcome would be very interesting since 
one would have found at least one example for a lattice regularized theory where fermions 
exhibit a chiral coupling to an "external" gauge field. In contrast in four dimensional models 
the naive power counting analysis suggests that the coupling to the effective gauge field U'^^ 
vanishes like a^ in the continuum limit. 



The U(l) pure gauge model ( |2.2[ ) is confining for all values of the gauge coupling g. The 
confinement phase is expected to be present also in the two dimensional U(l) Higgs model at 
small values of k and to turn in the limit (7 — > into the vortex phase. For g > the scalar 
particles get confined into massive bosonic particles. The effective gauge field in eq. ( [7.5|) 
can be written in the form 



K. 



z' + E^^ + iW, 



fix 



(7.6) 

where z'^ is given in eq. ([4.5|) and Hfj^^ and W^^ are interpolating fields for bosonic bound 
states in the confinement phase with quantum numbers J^'" = O"*"^ and 1 in lowest spin 
state ||3^. The field Hf^x couples primarily to the Higgs-like scalar particle according to 



H, 



/IX 



TTiH Hx 



(7.7) 



where mn is some mass scale which has to be introduced since the scalar field H^ is dimen- 
sionless (-ff^x is not a vector under lattice rotations, which forbids a relation of the form 
H^x -^ d^H{x)). The field Wfj,x couples to the vector boson. 

After inserting ( [7.6|) into ( |7.4] ) and a trivial rescaling of the fields \E')^ and \E'^ we obtain 
for Cp the form 



Cf 



■E 



^X l/J'^X+fl ^X+fM IfJ-^X 






mnl 



M=i 



z^ 2 
1 1 ^ 



7W-, .tM) vf7(") ^, iT/(") 



9 Z-^ ^-^x ^ L,x ifJ.^ L,x+/i ^ L,x+ix ifJ.^ L,x 



^ L,x ifJ.^ L,x+/i ' ^ L,x+/i ifJ.^ L,x 



(7.9) 
(7.10) 



The term ( [7.8|) describes a free neutral fermion with mass rrip = y/z. The expression ( [7.9|) 
suggests that the coupling of the neutral fermion to the Higgs-like bound state vanishes like 
a. However, the neutral fermion couples chirally in ( [7.10| ) to the vector boson field W^x if 
its dimension is one, as suggested by the naive dimensional analysis. 



In order to find out whether the fields W^x and H^x are indeed dimension one operators 
we have computed the scale dependence of the corresponding wave-function renormalization 



Figure 6: The wave- function renormalization constants Zw (squares) and Zh (circles) are 
plotted respectively as a function of m? = rriy^ and m = m\ for the fixed ratio mu/mw = 
1.14. The two dotted lines through the origin are drawn to guide the eye. 

constants Zh and Zw- From the naive dimensional analysis these wave-function renormal- 
ization constants Zh and Zw are expected to vanish like a^. To see whether this expectation 
is correct we have computed the momentum space propagators 

Gh{V2) = (^Y.Hi.H,ye'P'^^'-yA , Gw{p2) = (^ E ^1-^1.^*'''^"'"''^ (^-H) 

\ x,y I \^ x,y I 

in the confinement phase of the U(l) gauge-Higgs model for several values of g and k 
and fitted the results for sufficiently small p2 to the free boson propagator ansatz given 
in eq. (|5.5|), which for Gw is considered as a special case (pi = 0) of a free massive vector 
boson propagator, {5^y-\-p^py/m?'){rn?+p'^)~^. In fig. 6 we have displayed the resulting wave- 
function renormalization constant Zw (squares) and Zh (circles) respectively as a function 
of m? = o?rn^r^pf^ys and m? = a'^Tn?H,phys ^^ ^^e fixed ratio mH/mw = 1.14. In both cases 
the wave-function renormalization constants obey nicely a linear dependence, supporting our 
expectation from the dimensional analysis. 

This result suggests that the neutral fermion exhibits in two dimensions indeed a non- 
vanishing chiral coupling to the massive vector bosons. This coupling is universal in the 
quenched approximation since the field W^x is proportional to a current which is conserved 
in the two dimensional U(l) gauge-Higgs model. The properties of this coupling in the full 
model with dynamical fermions are, however, not yet clear to us. 



The fermion couplings in the VXS phase may then be summarized quahtatively by the 
following effective lagrangian 

C^f = V^(")#^"^ + mP^^"-^^^-^ + ^^?i"V,4"%(^) , (7.12) 

where g^ = \/ Z-\^ j z^ and W'^^ is the vector field in the continuum with standard normal- 
ization. This expression for the effective lagrangian gives a satisfactory description of the 
fermion couplings at distances which are large in comparison with the typical length scale 
of the vector boson bound state. When lowering the value of the gauge couplings g the 
string tension becomes smaller and the bound states extends over larger distances. At small 
distances the scalar particles are then almost free and a more appropriate form of the action 
is given then by 

Cf -. ^("V^(") + mi.")v^^"V^"^ + \^fl^i^P{<\>*d,<\> - <pd,<P*) . (7.13) 

This form is expected to describe the fermion couplings in particular in the global limit, 
i. e. g = 0, where the bosonic spectrum consists of unbound scalar particles. The massive 
neutral fermion in eq. (|7.13|) is coupled to a two particle current. This interaction seems, 
however, to be quite weak since our numerical results for the neutral fermion propagator are 
in good agreement with the analytic prediction from the lowest order hopping expansion, 
which leads to a free fermion propagator. Furthermore our numerical results for the binding 
energy e^ suggest that the interaction in (|7.13| ) is too weak as to give rise to the formation 
of a <|)-\[f("') bound state. Also the expressions ( [7. 121) and (|7.13| ) are certainly different from 



the target model in eq. (p. 5]) which we had originally in mind. 



8 Conclusion 

We started out our investigations from the lattice lagrangian given in eq. ( p.l|) in the hope 
to obtain in the continuum limit the target model in eq. ( |2.5| ). In the global limit of the 
model {g — > 0) the unwanted species doublers can be removed completely from the spectrum 
within the strong coupling vortex phase (VXS). The physics in this phase differs, however, 
substantially from the target action which we had originally in mind. In the previous chap- 
ter we have developed two different scenarios (A and B) for the effective theory in the VXS 
phase which were summarized by the continuum lagrangians ( |7.3| ) and (|7.12 ), ( [7.13 ). Our 



numerical results are in favor of scenario B: In this case the fermionic spectrum in VXS 
contains only a neutral fermion which has zero charge with respect to the U(1)l group and 
which in the global limit g ^ exhibits a left-handed coupling to a two particle current. 
This coupling, however, appears to be weak since the neutral fermion propagator data are 
in nice agreement with the results from the lowest order hopping expansion which implies a 
free fermion behavior. Furthermore this coupling seems to be too weak as to give rise to the 
formation of a $-\l/(") bound state. When the gauge coupling is turned on, we argue that the 
neutral fermion couples chirally to the massive vector boson state in the confinement phase. 
If this scenario is correct, it would have been the first time that a chirally coupled fermion 
has been detected on the lattice. This result is also different from the previous findings in the 
strong coupling symmetric phase (PMS) of the fermion-Higgs models with Wilson- Yukawa 
coupling in four dimensions where the coupling of the neutral fermion to the bosonic bounds 
states vanishes presumably as a power of the lattice spacing |]I3 |. 
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